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Abstract. - We investigate how the normal energy transport is realized in one-dimensional 
quantum systems using a quantum spin system. The direct investigation of local energy dis- 
tribution under thermal gradient is made using the quantum master equation, and the mixing 
properties and the convergence of the Green-Kubo formula are investigated when the number 
of spin increases. We find that the autocorrelation function in the Green-Kubo formula decays 
as ~ t~^'^ to a finite value which vanishes rapidly with the increase of the system size. As a 
result, the Green-Kubo formula converges to a finite value in the thermodynamic limit. These 
facts strongly support the realization of Fourier heat law in a quantum system. 



The Fourier heat law is one of the most important properties in the nonequilibrium thermo- 
dynamics. It states that the heat current per volume is proportional to the thermal gradient. 
Microscopic dynamical origin for the realization of Fourier heat law has been actively studied 
using many Hamiltonian systems jl], ^, . In the complete harmonic chain, no global ther- 
mal gradient appears, and local equilibrium is not realized, which are attributed to the lack 
of scattering between modes (|, In an isotropic d-dimensional classical Fermi-Pasta-Ulam 
(FPU) system which has the nonlinear potential term, the mixing property satisfies because 
the auto-correlation function of the energy current roughly decays as t~'^l'^ . However due to 
its slow relaxation of the current fluctuation, the Green-Kubo formula diverges in the one- 
and two-dimensional systems which causes an anomalous energy transport |^,|^. 

In low dimensional systems, most of the problems preventing a normal thermal conduction 
arises from a slow fluctuation of energy current, including the failure of mixing property. 
This situation is also the case in the quantum systems. Many integrable one-dimensional 
systems show the failure of mixing property. In the isotropic Heisenberg chain, the energy 
current operator commutes with the Hamiltonian, which trivially causes the failure of mixing 
property and the abnormal energy transport [||. The recent experiments confirmed a such 
ballistic heat transport in Sr2Cu03 [Q and CuGeOs |^ which are described by the isotropic 
Heisenberg chain. 

In this paper, we study how the normal thermal conduction in a one-dimensional quantum 
system is realized in microscopic point of view. In quantum systems, the dynamical origins of 
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normal thermal conduction and the related quantum effects are not enough understood y,^. 
In classical systems, chaos characterized by the sensitivity to the initial state can be induced 
by nonlinearities and can play a crucial role for the mixing property. However such chaos 
cannot be induced by a linear Schrodinger equation in quantum systems pO| . In spite of such 
discrepancies in dynamics, the Fourier heat law can be realized in quantum systems. We 
show the strong evidence of the normal thermal conduction in a one-dimensional quantum 
spin system by direct investigation of properties under thermal gradient and the Green-Kubo 
formula. We find that the autocorrelation function in the Green-Kubo formula approximately 
decays as t~^'^. Thus the mixing property is satisfied and the formula converges. 

We consider the tractable simple quantum system whose Hamiltonian is described as, 

Ni N N 

n = Y^ Jalal^, + D {olal^, a^a^i^,) + T ^ + if ^ a|, (1) 

1=1 i=i 1=1 

where cr" (a = x, y, z) is the a component of the Pauli matrix at the fth site. N is the number 
of spin, and A^i is taken as — 1 for the open boundary condition and N for the periodic 
boundary condition, respectively. The first term is the nearest neighbor Ising interaction, and 
the second term is the Dzyaloshinsky-Moriya (DM) interaction, where the anisotropic vector 
is taken as D = (0, Z?, 0) for the general form D ■ [at. x Uf+i). The third and fourth term are 
the Zeeman term of x and z direction, respectively. 

This Hamiltonian shows the variety of symmetries by controlling parameters (J, I?, F, 
In the case of periodic boundary condition, there exists the translational symmetry (T„ <-> 
On+i {£ — 1, • • • , N) . In the absence of the DM interaction, i.e, D = 0, the system has the 
reflection symmetry, (T„ <-> a^-n- In the case of D = and H = 0, the system Hamiltonian 
is diagonalized by using free fermions through the Jordan- Wigner transformation . When 
F = and H = 0, the system has the time reversal symmetry which yields the Kramers 
degeneracies. In many systems with symmetries, conserved quantities prevent a normal heat 
transport, so that the Green-Kubo formula diverges |€,12|. 



As well known, level spectrum shows the universal statistics according to the noninte- 
grability of system We investigate the cumulative spacing distribution I{S) defined by 
I{S) — duP{u) for the energy spacing distribution P{u). We show the distributions in 
Fig.l for [J,D,T,H) = (0.2,0.2,0.2,-0.2) and (0.2,0,0.2,0) with the open boundary con- 
dition iVi = — 1. The former system has no trivial symmetries, while the latter can be 
mapped into free-fermion. Points are numerical data by diagonalizing the Hamiltonian (|l|), 
and the lines are theoretical predictions, i.e., the Poisson distribution I{S) = 1 — exp(— S*) for 
integrable systems, and the Wigner distribution I{S) = 1 — exp(— 7r5'^/4) for nonintegrable 
systems, respectively. The figure for the parameters (0.2,0.2,0.2,-0.2) shows the agreement 
between the numerical data and the Wigner surmise. For small spacing S* ^ 1, I{S) behaves 
as I{S) oc 5'i-92±o.02^ which is very close to the Wigner distribution I{S) oc S"^. This universal 
feature of spectrum indicates a complex dynamics of this system. Thus throughout this paper, 
we take the parameters as J = = F = —H = 0.2 

In order to investigate properties of the energy transport under thermal gradient, the 
reservoirs of different temperatures are directly connected with the system at the both ends. 
Here we adopt the phonon reservoir with the Ohmic spectral density for which the master 
equation of the system is written as (see e.g., and references therein); 

^ = [H, p{t)] A {C^p{t) + C^At)) , (2) 
where the first term in the right-hand side corresponds to the pure quantum dynamics of the 
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Fig. 1 - Cumulative level spacing distribution for A'' = 12 with the open boundary condition. The 
asterisk and cross are for the cases of (J, D, T, H) — (0.2, 0.2, 0.2, —0.2) and (0.2, 0, 0.2, 0), respectively. 
For small S in the former case, I(S) is approximately fitted by I{S) oc 5'i-92±o.02^ 



system, and £l and £r express the dissipative effects of the inverse temperature /3l at the left 
reservoir and /3r, at the right reservoir, respectively. The parameter A is the coupling strength. 

The dissipative term Ca (a — L,R) is given by Cap{t) — {\^Xa, Rap{t)] + [Xa,Rap{ty\^^ , 
where and are the system's operators directly attached to the left and right reservoir, 
respectively. Here we take Xl = crl and Xr = (t|^. The operator is given by {k\Ra\'m) = 
(i?fe — £',„)(e'^"^^'=~'^™^ — 1)^'^ {k\Xa\m) in the representation diagonalizing the Hamiltonian 
(|l|) as 'H\k) = Ek\k) and ■H\m) = S„|m). 

We first consider the energy profile defining the ^th local Hamiltonian 7i^; 



J2 i^^k + Hal) 



(3) 



k=l 



We numerically integrate the equation (H) and obtain the stationary density matrix pst when 
the initial density matrix is taken as the canonical distribution of temperature 1//9r. The 
simulation was carried out by the fourth order Runge-Kutta method with the time step 0.01 
for the parameters /?l = 0.5, /?r = 0.2, and A ~ 0.01. In Fig. 2, we show the energy profiles at 
the stationary state Tr {pstHi) for iV = 7, 8, 9, and 10. The two large circles are equilibrium 
energy values of inverse temperatures /3l and /3r, respectively. iV = 10 was the largest number 
within the present computational ability. The figure shows the profiles with a finite gradient 
which is the sufficient condition for the Fourier heat law, although we cannot enumerate the 
energy profile and energy current in the thermodynamic limit. The local energies at the both 
edges are different from the expectation values of the temperatures of reservoirs. This is 
attributed to the small coupling strength A(= 0.01). Actually we confirmed that for larger 
A (<C J, 1?), the energy profile smoothly changes and these differences become smaller. We 
found the qualitatively same finite energy gradient in the case of smaller temperature difference 
= 0.3 and /3r = 0.2. 

We next study the local energy distribution at the stationary state focusing on whether 
the local equilibrium is realized or not for the system oi N = 10. The local equilibrium can 
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Fig. 2 - The energy profile at tlie stationary state for from A'' = 7 to TV = 10. The abscissa is the 
scaled spin number, l/N, I — 1, ■ ■ ■ , N — 1. 



be realized even in finite size systems [|15|,|6|. We study the energy distribution of subsystems 
composed of two spins and four spins. We consider four local subsystems with two spins, 
i.e., the subsystems described by the local Hamiltonians, 7^2, '^4, '^6i and Tig, and two local 
subsystems with four spins, i.e., the subsystems described by {Ti.2 + Hs) and {He + Hr). 
Temperatures at boundaries are set as (/3l, /9r) = (0.5, 0.2) for the former case and (/3l, /3r) = 
(0.3,0.2) for the latter case. The £th local energy distribution Pi{ei) is calculated using the 
£th local reduced density matrix pi as, 

Pei^i) = {ie\pi\ii), (4) 

where \ie) is the ith eigenstate of the ^th local subsystem Hamiltonian. The numbers of eigen- 
states are 4 and 16 for the cases of (/3lj/3r) = (0.5,0.2) and (/3l,/3r) = (0.3,0.2), respectively. 
The £th local density matrix pe is obtained by taking the trace for stationary density matrix 
Pst in the Hilbert space exclusive of spins which belong to the ith subsystem. In Fig. 3, the 
energy distributions Pi{ei) are shown on the semi-log scale. The points linked by one line are 
the distribution at one local subsystem. With the increase of the subsystem number £, the 
gradient of the line becomes smaller. The distributions approximately take the exponential 
form. Thus this monotonic £ dependence of gradient of energy distribution indicates that 
this system can be the candidate which satisfies the local equilibrium. We must note that 
the Ising and DM interaction between nearest neighbor local subsystems and the existence of 
finite energy flow can cause deviations of Pe,{£i) from the exponential form. 

Generally even if the energy profile is well scaled and the local equilibrium is realized, 
thermal conductivity can diverge in the thermodynamic limit due to the slow relaxation of 
energy currents mentioned in the introduction ||l^ . We now investigate whether this type of 
divergence occurs or not in the present system. The Green-Kubo formula is derived on the 
basis of the local equilibrium state. Since the present system is the candidate which satisfies 
the local equilibrium, we expect that the Green-Kubo formula quantitatively describes the 
thermal conductivity k(/3) in the thermodynamic limit. The Green-Kubo formula reads as, 

<f3) - du{[jj{u)}), (5) 
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Fig. 3 - Local energy distribution for size A*' = 10. (a) Four local subsystems composed of two spins 
for (/3l,/3r) = (0.5,0.2). (b) Two local subsystems composed of four spins for (/3l,/3r) = (0.3,0.2). 



where {., .} means the anti-commutation relation, and (...) is the equilibrium average at the 
inverse temperature (3. The operator J(u) is the total current operator at the time u in 
the Heisenberg picture, i.e., J{u) — exp{iTi.u)Jexp(—iTiu). The total current operator J is 



calculated by the continuity equation of energy as J 



L, hi 



he+2, he 



using 



D {afa^^. 



limit, we define A{t) and C{t) as follows. 



^l^i+i) H'^i- order to consider the thermodynamic 



A{t) = ^({j, J(t)}), C{t)=f^ duA{u), 



(6) 



where A{t) is the autocorrelation function for the total current, and C{t) is the integration of 
A{t). We call C{t) the Green-Kubo integral. The numerical calculations are carried out for 
P = 0.3 with the periodic boundary condition iVi = N . In Fig. 4(a), we plotted the absolute 
values of the autocorrelation functions A{t) for system sizes N — 12, 14 and 16 on the log- 
log scale. Up to TV = 14, we used the numerical diagonalization for the Hamiltonian, and 
for iV = 16 we calculated the time evolution of wave function for the randomly chosen 1024 
initial states |9|. The figure indicates the power law decay ~ t^^-^ of Ait). In the inset, the 
raw data of A{t) for various system sizes {N — 8, 10, 12, 14, and 16) are plotted on the normal 
scale. All data are saturated to some finite values with fluctuations. This constant is partially 
caused by the energy degeneracies due to translational symmetry of periodic condition. Thus 
the behavior of A{t) is roughly represented as follows. 



A{t) - a{N,t)t- 



-1.5 



b{N,t) + B{N), 



(7) 



where a{N, t) and h{N, t) are some fluctuating functions with a vanishing mean value and B{N) 
denotes the saturated average value for N . B{N) can be exactly calculated by a numerical 
diagonalization up to TV = 14. For N = 16, it is calculated by directly averaging the numerical 
data of A{t). We flnd that B{N) vanishes with the increase of N following or faster than the 
exponential function 0.1 x exp(— O.SA^), as shown in Fig. 4(b). There we plotted the data 
from A'' = 6 to = 16. Thus the integral C{t) is expected to converge showing roughly the 
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Fig. 4 - (a): The absolute values of autocorrelation function as function of time for A'^ = 12, 14, and 
16. In the inset, the raw data are shown for N — 8, 10, 12, 14, andl6. (b): The average constant 
B{N) as a function of A''. The dashed line is the exponential function 0.1 x exp(— 0.5 x A^). 



following behavior in the thermodynamic limit iV — s- cx), 

C{t)r^j a{N,u)u-^-''' du + BiN)t ^ J a{oo,u)u-^-^ du < Cq + r^'^ . (8) 

In Fig. 5, we present the numerically calculated C{t) for various N . The behavior of C(t) 
strongly supports the convergence of the Green-Kubo integral in the thermodynamic limit. 
The behavior of C{t) reminds us the order of limit, i.e., taking the limit N to infinity first 
and then t to infinity. From the convergence of C{t), we conclude that the present one- 
dimensional spin system without trivial symmetries shows the finite thermal conductivity 
which is independent of the system size as far as N is very large. 

In conclusion, wc find the strong evidence of normal thermal conduction in the present 
system by microscopic investigation. The autocorrelation function in the Grcen-Kubo formula 
shows power law decay t~^'^ , which may be the characteristics in the high temperature 
region, i.e. semiclassical region It is also interesting to study the behavior at very low 
temperatures. The Fig. 3 indicates that this system can be the candidate which satisfies the 
local equilibrium. We should carefully and systematically investigate how the exponential 
form in the local energy distribution are realized. 

In the recent experiments showing a ballistic heat transport in of Sr2Cu03 0| and CuGeOa 
1^ , the energy is transmitted by spin excitations (magnon), and the coherent length is very long 
(about 100 times lattice constant in CuGeOs above the spin-Peierls temperature). Theoretical 
studies for temperature dependence of thermal conductivity are actively done |jl^,|l^ . It is also 
interesting to investigate theoretically and experimentally how the coherent length changes 
when the symmetries disappears in these realistic systems. 
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Fig. 5 - The Green-Kubo integral C(t) as a function of time for N — 8, 10.12, 14, and 16. 
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